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Abstract. - The atomic structure of self-assembled quasi-one-dimensional Gd chains on Si(lll) 
has been investigated by low-energy electron diffraction and scanning tunneling microscopy. Based 
on comparison between Gd and Ca chains we show that this Gd induced surface reconstruction 
belongs to the class of honeycomb chain-channel structures. This clearly demonstrates that, 
besides monovalent and divalent adsorbates, also trivalent adsorbates such as Gd stabilize silicon 
honeycomb chains. Consequently silicon honeycomb chains emerge as an universal building block 
in adsorbate induced silicon surface reconstructions. 
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Introduction. — Self- assembled atomic chains on sil- 
icon surfaces have been the focus of intense research be- 
cause of their quasi-one-dimensional (ID) electronic prop- 
erties and their interesting physics. Recently the fluctua- 
tion and condensation phenomena at the metal insulator 
phase transition of the In/Si(lll) system could directly be 
visualized via scanning tunneling microscopy (STM) [1-4]. 
Competing periodicities in fractionally filled bands lead to 
the coexistence of different Peierls distortions for the gold 
induced reconstructions [5-7]. 

Another important class of ID systems are the alkali 
metal (AM=Li, Na, K, Rb, Cs) and Ag induced, in- 
sulating (3x1) reconstructions formed by the deposi- 
tion of 1/3 monolayer (ML) onto the Si(lll) surface. 
The AM/Si(lll) systems adopt the so-called honeycomb 
chain-channel (HCC) structure [8-10] shown in Fig. []Ji) 
which is stabilized by the transfer of one electron from the 
monovalent AM adsorbate into the Si surface states. 
A very similar reconstruction with (3 x " 2" ) periodicity is 
formed by adsorption of alkaline-earth metals (AEM=Mg, 
Ca, Sr, Ba), where the x"2" notation stands for a x2 pe- 
riodicity along the adsorbate chains but missing coherence 
between adjacent chains [11]. Due to the divalency of the 
adsorbate only 1/6 ML, i.e. half the AM coverage, is re- 
quired to stabilize the HCC structure [12]. At 1/2 ML 
divalent adsorbates induce a (2x1) phase which was pro- 
posed to be formed of 7r-bonded Seiwatz chains shown in 
Fig. [TJ)) [13,14]. For intermediate coverages, a series of 
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Fig. 1: (Color online) a) Honeycomb chain and b) Seiwatz 
chain model with adsorbates lying in the channels between the 
chains. To stabilize the honeycomb chains, monovalent atoms 
are required to occupy every site in the channels (blue and 
yellow circles), whereas divalent adsorbates only occupy every 
second site (blue circles only) . Seiwatz chains are stabilized by 
divalent adsorbates. Open circles are Si atoms. The (3x1) and 
(2x1) unit cells are also shown. 



ID (nx"2") reconstructions, with n taking the values 5, 7 
and even 9 depending on the adsorbate, is formed which 
are considered to be composed of an appropriate combina- 
tion of honeycomb chains and Seiwatz chains (see Fig. [3] 
for the 5x"2" case). Similar series of reconstructions were 
also observed for the divalent rare earth metals (REM) 
Sm, Eu and Yb [15]. These REMs more commonly occur 
in the 3+ valence state, but depending on their chemical 
surrounding the 2+ configuration is occasionally preferred 
as in this case. Thus up to now, only monovalent and diva- 
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lent adsorbates were found to stabilize Si reconstructions 
containing the honeycomb chain building block. 
In this letter we focus on trivalent REMs, which ex- 
hibit chain structures with (5x"2") periodicity only, 
but whose detailed atomic structure has not been in- 
vestigated. Combining low-energy electron diffraction 
(LEED), STM and recent angle-resolved photoemission 
spectroscopy (ARPES) results [16], we show for the first 
time that the structure induced by trivalent adsorbates 
contains the same honeycomb and Seiwatz chains as in 
the chain reconstructions induced by divalent adsorbates. 
The use of multiple complementary surface analysis tech- 
niques is mandatory in the present case in order to derive 
a reliable structural model. Based on electron counting we 
are also able to explain, why only the (5x"2") periodicity 
is stabilized for trivalent adsorbates. 

Experiment. — We choose to investigate the Gd sys- 
tem, since it has recently been demonstrated that predom- 
inantly single domain atomic Gd chains can be grown on 
stepped Si(lll) having a slight misscut of 1.1° towards 
the [112] direction [17] allowing the use of macroscopic 
diffraction methods without domain averaging. Qualita- 
tively similar results are expected for the observed (5 x " 2" ) 
reconstruction induced by other trivalent rare-earth met- 
als Dy [18], Er [19] and Ho [20]. Gd was evaporated from a 
water cooled e-beam evaporator with of flux of 0.5 x 10 -4 
ML/s at a pressure below 5 x 10 -10 mbar onto the clean 
Si(lll)-(7x7) substrate held at 680 °C. The substrate was 
heated by passing a direct current along the step direction 
[110]. Growth and experiments were carried out in an ul- 
tra high vacuum chamber with a residual gas pressure of 
3 x 10 -11 mbar equipped with an Omicron LT-STM and 
Omicron Spectaleed LEED/ Auger optics. For STM mea- 
surements we used etched W tips. 

Results and Discussion. — Figure [2^) shows the 
LEED pattern of a typical Si(lll)-(5x"2")-Gd surface 
with one dominant domain and insignificant contributions 
from the two others and the Si(lll)-(7x7) reconstruction. 
Only the (5x1) spots sketched in Fig. (2}d) are clearly vis- 
ible. The x2 periodicity along the chains manifests itself 
through faint half-order streaks parallel to the x5 spots 
(not shown) observed only at certain energies. Similar 
streaks were reported in studies of divalent adsorbate sys- 
tems and explained in terms of a stochastic distribution 
of adjacent chains with random registry shifts leading to 
a (5x"2") spot pattern with its characteristic weak half- 
order streaks [11,21,22]. 

Whereas the LEED spot positions only determine the type 
of Bravais lattice of the surface structure, i.e. its trans- 
lational symmetry properties, the point symmetries can 
be determined by a symmetry analysis of the intensity vs 
voltage (IV) curves. The threefold symmetry of the unre- 
constructed Si(lll) surface termination is broken by the 
growth of the chains. Whereas the (0,-1) and (-1,1) beams 
are still equivalent as for the substrate, the (1,0) beam ex- 
hibits a distinctive spectral signature as can be seen from 
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Fig. 2: a) LEED pattern of Si(lll)-(5x"2")-Gd at 48 eV. b) 
Sketch of the 5x1 LEED pattern in reciprocal space with beam 
indices. In real space the chains are running along the vertical 
axis, c) Comparison between experimental LEED IV curves of 
Si(lll)-(5x"2")-Gd (full lines) and Si(lll)-(5x"2")-Ca (dot- 
ted lines). 



Fig. [2k). Thus only a mirror plane perpendicular to the 
chains is retained. 

To obtain information about the atomic positions we com- 
pare LEED IV curves from Si(lll)-(5x"2")-Gd to the 
curves from Si(lll)-(5x"2")-Ca in Fig. Et). IV-LEED 
fingerprinting has played a crucial role in establishing the 
equivalence between different AM induced (3x1) HCC re- 
constructions, since it was recognized that the Si(lll)- 
(3xl)-AM reconstruction is predominantly a substrate re- 
construction with a common structure independent of the 
adsorbate species [23]. Visual inspection of Fig. [2t) al- 
ready shows that the agreement between the Gd induced 
and the Ca induced reconstruction containing one honey- 
comb chain and one Seiwatz chain is surprisingly good. 
Most peak positions of the Gd chains fall on the same en- 
ergies as for the Ca chains with comparable relative inten- 
sities. To obtain a quantitative measure for the agreement 
between the two structures we calculated Pendry's R fac- 
tor R p [24], which takes into account the peak positions 
but also the relative intensities between the peaks. For 
the integral order spots we obtain R p = 0.29. For the 
fractional order beams we obtain R n = 0.35. These val- 
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Fig. 3: (Color online) Structural model for the Si(lll)- 
(5x"2")-Gd surface. Open circles are Si atoms, filled blue 
circles are Gd atoms. The (5x2) unit cell divided into two 
(3x1) and two (2x1) unit cells is also shown. Arrows indicate 
a registry shift of the adsorbates in the channel. 



ues are similar to R p = 0.36 obtained by Lottermoser [8] 
comparing theoretical curves to experimental data for the 
HCC model. The good agreement between the two ex- 
periments suggests that both structures share the same 
structural building blocks. Deviations may be due to the 
difference in atomic number and the associated scattering 
characteristics between Gd and Ca, differences in the pre- 
cise adsorption geometry and coverage. 
Adsorbate coverage is an important parameter for the de- 
termination of any structural model. The exact amount 
of Gd at the surface is difficult to determine accurately 
due to the fact that Gd diffuses into the bulk above 600 
°C [17]. The ideal adsorbate coverage can however be de- 
termined when considering the electron count required to 
stabilize the honeycomb and Seiwatz chains. The HCC 
structure is known to be stabilized by the donation of one 
electron per (3x1) unit cell [10,12]. Similarly the Seiwatz 
chain requires two electrons per (2x1) cell, since it may 
be stabilized by 1/2 ML of divalent adsorbates. This is 
consistent with the number of surface states observed in 
ARPES [25]. The (5x2) unit cell can be thought of as 
being build from two (3x1) and two (2x1) cells, thus re- 
quires six electrons to be stabilized. Since Gd is trivalent, 
the ideal coverage is two Gd atoms per (5x2) cell or 1/5 
ML. This is in agreement with the estimate of 0.2-0.4 ML 
given in Ref. [17]. 

The proposed model for the Si(lll)-(5x"2")-Gd surface is 
shown in Fig. [3] consisting of alternating hexagonal hon- 
eycomb chains and zig-zag Seiwatz chains made of Si. The 



adsorbates are expected to form chains in the channels in 
between. Due to the weak sensitivity of IV-LEED to the 
adsorbate itself, we can not decide which absorption site is 
favored. Any structural model must be consistent with re- 
sults from other experimental techniques. Fig. [4] presents 
STM images of the Gd chains. The overview a) shows 
long, parallel chains running along the [110] direction. The 
separation between the rows is consistent with the x 5 pe- 
riodicity observed in LEED patterns. High magnification 
empty and filled state images acquired in the same scan to 
preserve their mutual registry are shown in Fig. [4jo) and 
c) respectively. The structural model is superimposed. 
Based on simulated STM images derived from local den- 
sity approximation (LDA) calculations for the HCC struc- 
ture [10, 12], we identify the dark rows in the empty state 
image with the location of the honeycomb chains. High 
intensity in the empty state image is found along the ad- 
sorbate channels for both the honeycomb and the Seiwatz 
chain structure [12,26]. This is easily understood by notic- 
ing that the empty orbitals are necessarily located on the 
adsorbate atom, since it donates its electrons to the sili- 
con chains. The filled state image c) appear as triple rows 
of protrusions with x2 periodicity along the rows. The 
third row located along the Seiwatz chain (marked by S 
in Fig. [4k) appears to lie slightly lower than the two main 
rows (marked by H in Fig. Hh), which we identify with 
the honeycomb chains. In a previous STM study only the 
two main rows H were resolved [17]. The pairing of pro- 
trusions causing the x2 periodicity along the chains has 
been found to be rather electronic in origin than geomet- 
ric [27]. The electrostatic attraction between a positive 
adsorbate ion and the electrons in the neighboring satu- 
rated dangling bonds give rise to such paired protrusions. 
We also remark that the registry of neighboring chains is 
correct in our model. Careful inspection of the filled state 
STM image shows that the honeycomb chain comes in two 
configurations, either as two parallel rows of protrusions 
or in a zig-zag configuration, indicated by empty circles 
in Fig. [It). Such a registry shift of only one period be- 
tween the two rows of the honeycomb chain is illustrated 
by the arrows in Fig. [3] and is simply due to a missing ad- 
sorbate and a consecutive shift of all the following adsor- 
bates by one period along the chain direction. The local 
mixing of these two arrangements with poor long range 
order is responsible for the x2 streaks seen in LEED pat- 
terns [21]. Furthermore this kind of defect leads to a local 
charge imbalance. It has been suggested that additional 
Si adatoms are able to supply electrons that dope the par- 
ent chain structure [28] and may be able to compensate 
for such missing charge. Additional Si atoms are necessar- 
ily present since the formation of the HCC structure and 
consequently also of the (5x"2") structure is accompanied 
by significant Si mass transport at the surface [29] due to 
the fact that the Si atom surface density of the (5x"2") 
structure is not equal to that of Si(lll)-(7x7). Although 
steps may serve as a reservoir for reintegrating ejected Si 
atoms into the surface, electromigration due to dc current 
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Fig. 4: (Color online) a) STM topography overview (U= 1.9 V), 
60 nm x 90 nm, b) and c) high resolution topography of empty 
(U=1.9 V) and filled states (U=-1.9V), 18 nm x 9 nm, 1=0.18 
nA. Arrows indicate the location of the honeycomb chain (H) 
and Seiwatz chain (S). Empty circles mark the two possible 
configurations of the honeycomb chain caused by a registry 
shift of the Gd atoms in the adjacent channel as indicated by 
arrows in Fig. [3] 



heating parallel to the steps does not favor the Si atoms 
to wander towards the steps, resulting in a large number 
of randomly distributed protrusions on top of the chains. 

We now turn to the discussion of recent ARPES results 
from Si(lll)-(5x"2")-Gd [16], which provide additional 
confirmation for our structural model. At least three semi- 
conducting surface states are observed at binding energies 
between 1 and 2 eV, whose dispersions, band widths and 
symmetry properties are very similar to those of the AM 
and AEM induced (3x1) and (3x"2") reconstructions [30] 
supporting a honeycomb chain based structure. Further- 
more ARPES data for the AEM induced (5x"2") struc- 
ture resembles the one from the (3x"2") reconstruction. 
Very weak intensity is observed at the Fermi energy, but 
has been interpreted as being due to defect states. A small 
contribution to the spectral weight at the Fermi energy 
was also observed in the semiconducting Si(lll)-(3x"2")- 
Ca system [21], but was attributed to remaining (7x7) 
regions of pure silicon. Additionally, prolonged anneal- 
ing of the Gd induced reconstruction at 680 °C leads to 
the nucleation of metallic Gd silicide islands at the ex- 
pense of the chain reconstruction, which might possibly 
be responsible for the observed photoelectron signal at the 
Fermi energy. However, from STM measurements we do 
not find evidence for a metallic surface state localized on 
the chains. We conclude that the Gd induced structure is 
semiconducting and consequently requires an even num- 
ber of valence electrons per unit cell in agreement with the 
coverage of two Gd atoms per (5x2) unit cell. Therefore 
all ARPES results fully support our structural model. 
A peculiar experimental finding to be explained is that Gd 
and other trivalent REMs stabilize chain structures with 
the (5x"2") symmetry exclusively, whereas the monova- 



lent adsorbates stabilize only the genuine (3x 1) HCC 
structure and the divalent adsorbates induce a series of 
(nx"2") reconstructions. Monovalent adsorbates must oc- 
cupy every site along the channel between the honeycomb 
chains to satisfy the doping criterion. For lower coverages 
only parts of the Si(lll)-(7x7) are transformed, whereas 
higher coverages induce different surface structures. The 
stabilization of Seiwatz chains requiring two electrons per 
unit cell is not possible. Divalent adsorbates in turn must 
occupy every second site to satisfy the doping balance. 
For higher coverages however, additional adsorbates may 
be incorporated in the channels at the expense of reducing 
every second honeycomb chain into a Seiwatz chain. For 
trivalent adsorbates, charge balance requires that every 
third site in the channel is occupied, if one wants to build 
a structure exclusively formed by honeycomb chains. This 
is apparently energetically unfavorable compared to an oc- 
cupation of every second site, which requires the combina- 
tion of a honeycomb chain with a Seiwatz chain resulting 
in the (5x"2") symmetry. The stabilization of a (5x"2") 
period requires a total of six electrons, a condition easily 
satisfied by taking two trivalent adsorbates per unit cell. 
(7x"2") and (9x"2") reconstructions are not observed for 
the trivalent adsorbates. Consisting of one honeycomb 
chain and two respectively three Seiwatz chains, they re- 
quire 10 respectively 14 electrons per unit to be stabilized, 
a condition which can not be satisfied by trivalent donors. 
Electron counting thus provides a simple intuitive picture 
for the occurrence of the various phases. 

Conclusion. — Driven by the elimination of dangling 
bonds and relief of surface stress, silicon surfaces recon- 
struct in strikingly diverse ways. Among the large vari- 
ety of adsorbate induced reconstructions, the honeycomb 
chain emerges as a most stable building block allowing 
maximum reduction of the surface energy. The fact that 
only silicon atoms participate in the formation of the hon- 
eycomb chains allows a variety of adsorbates to adopt the 
HCC structure by donating the correct number of elec- 
trons to the substrate. Combining the complementary 
strength of IV-LEED fingerprinting, STM and ARPES, we 
demonstrated for the first time that next to monovalent 
and divalent adsorbates, trivalent adsorbates are also able 
to stabilize the honeycomb chains. Based on a intuitive 
electron counting model, we are further able to explain, 
why only the (5x"2") symmetry is stabilized by trivalent 
adsorbates. Our conclusions allow to enlarge the range 
of honeycomb chain stabilizing adsorbates to the trivalent 
elements. 
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